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(54) Tide: EPILEPTIC SEIZURE PREDICTION BY NONLINEAR METHODS 
(57) Abstract 



This invention relates to methods and apparatus for auto- 
matically predicting epileptic seizures and to monitor and analyze 
brain wave (EEG or MEG) signals. Steps include acquiring the 
brain wave data from the patient (11); digitizing the data; obtain- 
ing nonlinear measures of the data via chaotic time series analy- 
sis tools; obtaining time serial trends in the nonlinear measures; 
comparing the trends to known seizure predictors; and providing 
notification that a seizure is forthcoming. 
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EPILEPTIC SEIZURE PREDICTION BY NON-LINEAR METHODS 



The United States Government has rights in this invention pursuant to contract no. DE- 
AC05-84OR21400 between the United States Department of Energy and Lockheed Martin 
Energy Systems, Inc. 

FIELD OF THE INVENTION 

The present invention relates to the application of chaotic time series analysis (CTSA) to 
electroencephalogram (EEG) data and magnetoencephalogram (MEG) data, and more particularly 
to the analysis of the data to predict epileptic seizures. 

CROSS-REFERENCE TO RELATED APPLICATIONS 
This invention is related to Docket No. 1742-X Epfleptic Seizure Detection by Non-linear 
Methods by Lee M. Hively, Ned E. Clapp, C. Stuart Daw, and William F. Lawkins and to Docket 
No. 1744-X Method and Apparatus for Extraction of Low-Frequency Artifacts from Brain 
Waves for Alertness Detection by Ned E. Clapp and Lee M. Hively, both of which are filed on 
even date herewith, and both of which are assigned to the same entity. 

BACKGROUND OF THE INVENTION 
Hie theory of nonlinear dynamics provides a basis for understanding and poten tiall y 
controlling many complex physical and engineering systems. An extensive literature exists for 
nonlinear dynamics in the brain and related work (18). It is well known that brain waves exhibit 
seemingly random, unpredictable behavior, that is characteristic of deterministic chaos (1, 20). 
Moreover, chaotic behavior is "normal," while nonchaotic or periodic behavior is indicative of 
pathophysiology in experimental epilepsy (20). Schiff et aL (11) showed that chemically-induced 
seizures in rat-brain can be electrically controlled, leading to speculation (4) that human epilepsy 
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may be controlled without drug or surgical intervention. However, effective use of chaos control 
for epilepsy requires definitive seizure prediction. Thus, this invention diagnoses brain wave data 
via chaotic time series analysis (CTSA) methods to predict an epileptic seizure. 

Nonlinear analysis of neurological diseases via EEG oatt is extensive. For example, see 
the 1994 review by Elbert et aL (18). Epilepsy can be recognized only with clear EEG 
manifestations, but even these seizures are not easy to detect because there is no stereotyped 
pattern charaaeristic of all seizures (5). Work by Olsen and colleagues (7) used various linear 
measures with autoregressive modeling, discriminant analysis, clustering, and artificial neural 
networks. Valuable nonlinear tools for studying EEG data include correlation dimension, mutual 
information function, Kolmogorov entropy, phase-space attractors, and largest Lyapunov 
exponent. No publications or patents apply one or more CTSA measures to EEG data for 
systematic characterization of non-seizure, seizure, and transhion-to-seizure, for prediction of an 
epileptic seizure. 

Very recent analysis by TheileT (16) studied correlation dimension and Lyapunov 
15 exponent, using a form of surrogate analysis on a single EEG time series during an epileptic 
seizure. The surrogate analysis involved a random shuffling of blocks of time serial data, each 
block containing one quasi-periodic spike-wave complex. The auto-correlation function for the 
lg inal data is nearly ^distinguishable from the surrogate data. The correlation dimension for the 
,gmal data is significantly different from the surrogate data only at large scale sizes and large 
embedding dimensions. The maximum Lyapunov exponent (X) was negative for both the original 
and surrogate data and not substantially different, contrary to previous work which found positive 
X values. Theiler concluded that his analysis suggests a nonlinear oscillator with noise on the time 
scale of the spike-wave complex, but cannot indicate whether chaos exits on a shorter time scale. 
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OBJECTS OF THE INVENTION 

Accordingry.it is«iobje«ofthe present invention to provide methods and apparatus for 
predicting epileptic seizures in a patient and providing notification to the patient or to persons 
who can assist the patient, or to a data collection or transmission means. 

Further and other objects of the present invention will become apparent from the 
description contained herein. 
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SUMMARY OF THE INVENTION 

In accordance with one aspect of the present invention, the foregoing and other objects 
are achieved by a method for automatically predicting an epileptic seizure in a patient which 
comprises the steps of: providing at least one channel of a patient's raw brain wave data, called e- 
data, selected from the group consisting of electroencephalogram data and 
magnetoencephalogram data; separating the e-data into artifact data, called f-data, and artifact- 
free data, called g-data, while preventing phase distortions in the data; processing g-data through 
a low-pass fiber to produce a low-pass-filtered version of g-data, called h-data; applying at least 
one measure selected from the group of consisting of the linear statistical measures minimum and 
maximum, standard deviation, absolute minimum deviation, skewedness, and kurtosis, and the 
nonlinear measures time steps per cycle, Kohnogorov entropy, first inmimum in the mutual 
information function, and correlation dimension to at least one type of data selected from the 
group consisting of e-data, f-data, g-data, and h-data to provide at least one time serial sequence 
of nonlinear measures, from which at least one indicative trend selected from the group consisting 
of abrupt increases, abrupt decreases, peaks, valleys, and combinations thereof is determined; 
comparing at least one indicative trend with at least one known seizure predictor; and 
determining from said comparison whether an epileptic seizure is oncoming in the patient. 

In accordance with a second aspect of the present invention, the foregoing and other 
objects are achieved by apparatus for automatically predicting an epileptic seizure in a patient 
which comprises: data provision means for providing at least one channel of raw brain wave data, 
called e-data, selected from the group consisting of electroencephalogram data and 
magnetoencephalogram data; separation means for separating e-data into artifact data, called f- 
data, and artifact-free data, called g-data, while preventing phase distortions in the data, 
connnunicabry connected to said data provision means; low-pass filter means for filtering g-data 
to produce a low-pass filtered version of g-data, called h-data, communicabry connected to said 
separation means; application means for applying at least one measure selected from the group 
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BRIEF DESCRIFnON OF THE DRAWINGS 
In the drawings: 

«. Jis obtained tan *« p*ta, digbixed, processed. and ansiyzed by nominee, -oefcods to 

predict epileptic seizures. ^ ^ _ y 

Rg. 2 sandatd EEG eteom.de positions on . patient's scalp for she topotar 

^auge, looking tan above, « odtad m Examples 1 and H. 

F , 3 snows satnple *~ of EEG data, to »«. saw data <e-d«. astbe dooed «rve 

„ eg * - — * * «■ - - — - — ^ " d 

^.ftbiskrventa. d——-- 
fig. 4 comprises Figs. 4.1, 4.2, 4..J, ana 

-hows linear and nonlinear measures of raw (e-data); 4.2 
ofartifact(f-data);4.3 shows linear and nonlinear measures 
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time serial data for Example L 4.1 &ows linear and nonlinear measures 
shows linear and nonlinear measures i 
ofarti&ct-filteTed(g-data);and4.4sl 

^(h-data). Various mathematical properties and character!*- are shown for each 
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data type, and are further shown in pans (a), (b), (c), (d), (e), (f), (g), (h), and (i) for each data 
type. These measures are shown as curves, each curve representing a time serial sequence of 
linear or nonlinear measures. Within each curve, significant features such as abrupt increases, 
abrupt decreases, peaks, and valleys may be viewed as indicative trends which are then compared 
to trends which have been shown to be siezure predictors. Thus from the comparison it can be 
determined whether a siezure is oncoming in the patient In parts (a), (b), (c), (d), and (e) for 
each data type, the solid line is the specific measure, the dashed line is the 1 1-point average, and 
the dotted line is the standard deviation for the measure. 

fig 5 comprises figs. 5. 1, 5.2, 5.3, and 5.4, and shows linear and nonlinear measures of 
time serial data for Example K 5. 1 shows linear and nonlinear measures of raw (e-data); 5.2 
shows linear and nonlinear measures of artifact (f-data); 5.3 shows linear and nonlinear measures 
of aitifart-fihered (g-data); and 5.4 shows linear and nonlinear measures of artifact- and low-pass 
filtered (h-data). Various mathematical properties and characteristics as computed for each data 
type are further shown in parts (a), (b), (c), (d), (e), (f), (g), (h), and (i) for each data type. These 
measures are shown as curves, each curve representing a time serial sequence of linear or 
nonlinear measures. Within each curve, significant features such as abrupt increases, abrupt 
decreases, peaks, and valleys may be viewed as indicative trends which are then compared to 
trends which have been shown to be siezure predictors. Thus from the comparison it can be 
determined whether a siezure is oncoming in the patient. In parts (a), (b), (c), (d), and (e) for 
each data type, the solid line is the specific measure, the dashed line is the 1 1-pomt average, and 
the dotted line is the standard deviation for the measure. 

For a better understanding of the present invention, together with other and further 
objects, advantages and capabilities thereof, reference is made to the following disclosure and 
appended claims in connection with the above-described drawings. 

25 DETAILED DESCRIPTION OF THE INVENTION 

Chaotic time series analysis (CTSA) is applied to human electroencephalogram (EEG) 
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and magnetoencephalogram (MEG) data. For reference, equation numbers are shown at the right 
of each equation. 

In Examples I and H, EEG data from three epochs were examined: epileptic seizure, non- 
seizure, and transition from non-seizure to seizure. The CTSA tools were applied to four forms 
of these data: raw EEG data (e-data), artifact data (f-data) via application of a quadratic zero- 
phase filter of the raw data, artifect-fihered data (g-data) that was the residual after subtracting f- 
data from e-data, and a low-pass-fihered version of g-data (h-data). The transition from normal 
to seizure state also is characterized by distinctly different trends in the nonlinear and linear 
measures, as seizure predictors. Several linear and nonlinear measures distinctly predkt an 
epileptic seizure, including sample standard deviation of the correlation dimension for e-data 
haviag small variation during non-seizure but undergoing large rises and fells during transmon 
from non-seizure to seizure; sample standard deviation for the correlation dimension of f-data 
having small variation during non-seizure but undergoing large rises and fells during transmon 
fromnon-seizureto seizure; sample standard deviation for the correlation dimension of g-data 
having small variation during non-seizure but undergoing large rises and fells during transmon 
fromnon-seizure to seizure; sample standard deviation of the correlation dimension for h-data 
having small variation during non-seizure but undergoing large rises and fells during transmon 
fromnon-seizure to seizure; the Kolmogorov entropy for f-data having large, aperiodic vacuous 
during transition from non-seizure to seizure and less variation during non-seizure; the 
Kolmogorov entropy for h-datahaving large, aperiomc variations durmg transmon from non- 
seizure to seizure and less vnrfetion durmg non-seizure; a gradual decrease 
for the Mutual Information Function for the e-data; a gradual decrease in the first mnumum for 
the Mutual Information Function in the g-data; the average Kolmogorov entropy for e-data nsmg 
m onotonicaily while the Kolmogorov entropy undergoes oscillations of increasing amphtude 
about the average |^^ta^^-'^ t ***■■ l,< * , 
variability in skewedness and kurtosis for f-data being moderate with small sample standard 
deviationsthenincreasing markedly with ^^^^nssa^^^ 
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skewedness returning to moderate values wide the sample standard deviation of the kurtosis 
remaining at 2 to 3 times the non-seizure value. Analysis of e-data shows that statistically 
significant nonlinear structure is present during the non-seizure, transition, and seizure epochs. 

Two sets of channel- 13 EEG data from one patient are provided as examples. In fig. 1 of 
the drawings, 11 shows the patient's head, looking from above. 13 shows EEG electrode 
positions on the patient's scalp. 15 shows a nonlinear measure of EEG data. m fig. 2 of the 
drawings, C13 labels the position where the channel 13 data, which is used in this work, 
originates. Analysis of C13 data was chosen for these examples to demonstrate the robust 
removal of eye-blink artifact, which otherwise dominates channel 13 because of its proximity to 
the eye. This method and apparatus can be applied to data from other EEG channels, as well as 
to MEG data, as is apparent to those skilled in the art Both sets of data included non-seizure, 
transition-to-seizure, and epileptic seizure data. The analysis included various linear measures 
(standard deviation, absolute average deviation, skewedness, kurtosis), phis nonlinear measures 
(time steps per cycle, Kolmogorov entropy, first minimum in the mutual information function, and 
15 correlation dimension). Four forms of the data were analyzed: raw EEG (e data), artifact (f data) 
via application of a zero-phase quadratic filter, artifact-filtered (g data) that was the residual after 
subtracting f-data from g-data, and a low-pass-filtered version of the g-data (h data). Both linear 
and nonlinear measures provide trends to predict seizure onset. Nonlinear tools identified 
predictive trends with and without artifact removal, with and without low-pass filtering, 
demonstrating the robustness of these methods to noise and artifacts. Surrogate analysis of e-data 
showed that this data has significant nonlinear structure. 

Analysis in accordance with the present invention reveals trends in linear and nonlinear 
measures that precede the seizure by 500-800 seconds, as shown in Figs. 4. 1-5.4, and as 
described further herein. These trends are easily computable and permit seizure prediction. 
25 Transition indicators are somewhat different for the two seizure examples. This indicates that 

aezure indicators that are chosen for comparison with indicative trends for siezure prediction in a 
patient must be chosen carefully. 
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channel, the brain-wave signal (e) at time (t) is sampled at regular intervals ft =i At) to yield 
of time serial data e, = eft). A fiher-window length was chosen of 2n + I points from the time 
series, where n is the number of points on cither side of the central point (ej as indicated in the 
sequence below. 

central point 

I 

e Mi» e e**i» — » Ce-n e e > e^,, ... , e^,, e^ 



n points n points 



The data was fitted to a quadratic equation that takes the form F ; = Fft) = a, ft - tj 2 + a, 
ft - 1 J + a, = a, T 2 + a 2 T ; + a,. Here, t e = c At is the time at the central point, and T; = t, - 1^ 
This approximation is fitted to the data, by minimizing the sum of squares of the differences 
between the quadratic equation, F(t), and the raw EEG data, eft), corresponding to the 
in the following function: 

1 = £ [*(',)-e(',)] 2 ■ E Kofi+frt-*^ Eq. 1 



The mmimumin Lis found from the condition 3U3a k = 0, for k ={ 1, 2, 3), formingthree 
simultaneous linear equations in three unknowns. The window-averaged artifact <FJ is given by 
the fitted value of the central point, F e = F(0) = a,. Note that the sums over odd powers of % are 
zero and that symmetric sums over even powers of T, (over i from -n to +n) can be convened to 
sums from 1 to n with T f = t At, yielding a window-averaged solution for the artifact signal: 

c — — : Eq. 2 

(4n 2 +4*-3X2/i+l) 

Here, 2 f indicates the sum over i from -n to +n. Sums over even powers of T were 
explicitly evaluated with standard formulae (10). The effort to evaluate F e can be reduced 
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substantially by computing the sums initially from Eq.2 (at c = n * 1), and then using the 
following recursions thereafter. 

n 



o 



Eq.3 



Eq.4 



Eq. 5 
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low-frequency component while retaining the higher frequency information. 

For a specific EEG channel, a time history of the nonlinear measures was obtained by 
applying the CTSA tools to a series of 20-second analysis-windows of the four data types (e, f 
g,h). These data are designated herein as x. The length of the analysis window (w) was 10,240 
points. Each analysis-window had a 5, 120-point overlap with the previous (or next) analysis 
window of data. This 50% overlap provides an optimal mix of new and old data for smooth time- 
history trend generation (2), as illustrated below. 



|< — o to 20 seconds — >| 

|<— -10 to 30 seconds — >| 

|< — 20 to 40 seconds — >| 



|<-100 to 120 seconds->| 

|<~110 to 130 seconds->| 



The zero-phase quadratic filter provides artifect-fihered data with frequencies of* 2 Hz. 
A heuristic for (linear) Fourier analysis is that * 10 periods of data are required to taithfuDy 
recover cyclic information at a specific frequency. Thus, *5 seconds of data are needed to obtain 
Fourier amplitude and phase information at a signal frequency of 2 Hz. However, this heuristic 
does not apply to nonlinear analysis. For example, 20 seconds of data are necessary to obtain 
consistent results for the Kolmogorov entropy. This need for longer dataset lengths (-10,000 
points) for consistent nonlinear measures conflicts with the need for shorter dataset lengths 
(s5,000 points) to provide adequate resolution for the time history generation of trends. 
Consequently, a 20-second analysis window was used, as described above, with a 50% overlap 
25 for an effective time history resolution of 10 seconds. The nonlinear measures for each 20-second 
analysis window were associated with the rime at the center of the analysis window, Le., every ten 
seconds. The nonlinear measures for each 20-second analysis window were associated with the 
time at the center of the analysis window, Le., every ten seconds. Shorter or longer analysis 
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window lengths can be used in proportion to higher or lower data sampling rates, as is obvious to 
those skilled in the art. 

Many charaaerization tools exist for chaotic data analysis. A subset of the tools that were 
found in previous work to be good measures for EEC data were used in carrying out the present 
invention These tools include the following: standard statistical measures (ixnnhnum, xnaxunum, 
average absolute average deviation, standard deviation, skewedness, kurtosis, time per cycle); 
Kolmogorov entropy and entropy spectrum; mutual information function; maximum KkeKhood 
correlation dimension and correlation dimension spectrum; surrogate generation and nonhneanry 
tests- and nonlinear digital filters (as discussed herein). 

' ^opy, correlation dimension, and mutual information were used in the present invention 

asnonlbearmeasuresfor seizure analyst The 

the time scale for generating the return map for EEG dynamics. The return map underlies the 
correlation dimension (measure of dynamic complexity) and entropy (measure of dynanuc 
predictabilhy). Wehave appliedthese same measures successfully in analyzing other system, 

The statistical measures for the present study are obtained by standard methods (6). The 
n ^ a ndmhumumareobtame4asn— and minimum (respectively) over the x, values 
in , time-serial window of w points. The average (x) is given by. 



The r-th order moment (m,) of toe x-datt is: 

h t. E"!' 7 

Toe ebsohrt« sver^e deviedoo (.) provides a roooa Wico.o, of roe x, variety sod is 
defined as: 

14 



WO 97/34525 



PCT/US97/04475 



a = 



f j \ 



Eq.8 



An unbiased estimate of the standard deviation (o) is: 



o = 



J/2 



w-1 



Eq.9 



An estimate for the skewedn ess (s) is: 



s - 



m 



3/2 



Eq. 10 



An estimate for the kurtosis (k) is: 



m 4 

k = ~r> ~ 3 Eq i1 



The average cycle time (T c ) is important as a characteristic time of the nonlinear system 

7. _ window length in timesteps 

' number of mean crossings] Eq. 12 



( 



■ 



The mutual information function (MIF) is a nonlinear version of the (linear) auto- 
correlation and cross-correlation functions, and was developed by Fraser and Swmney (3). 
Mutual information measures the certainty with which a measurement can be predicted, given the 
outcome of another related measurement. Examples of the later include the same EEG channel at 
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a different time, and another EEG channel at the same (or different) time. The MIF indicates the 
average information (in bits) that can be inferred from one measurement about a second 
measurement, and is a function of the time delay (number of time steps) between the 
cneasurements. Tne mutual informauon function also measures the nonlinear time dependent 
correlation in the same signal For EEG data, information decay in an individual channel 
(univariate MIF) indicates local time scale, as the average time lag ( tj - 1 ; ) that makes tftf 
independent of ^ and corresponds to the first 0 ocal)n 1 htinmm(M l ktin^cpsJmtheMff 
(3) For use herein, a minimum is defined as two successive decreases in the signal value, 

yield less consistent results. Tne MIF, I(Q,S), and system entropy (H) for two measurements (Q 
and S) are defined by. 

I{QJS) = I{S,Q) = H{Q)+H{S)-H(S,Q) 

H(S) = -5>s(s,)log[W] Eq ' 14 

i 

h(s,0 = -E^^NV"^ 1 Eql5 
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(14) with modifications for noise (13): 
/ \ ( 



Eq. 16 



where M is the number of randomly sampled point pairs, r^ is the normalized maximum-norm 
distance between the (randomly chosen) i-j point pairs as defined in Eq. 17 (below), and r„ is the 
normalized distance (scale length) associated with noise as measured from the time serial data. 
The distances are normalized with respect to some nominal scale length (Lo), Le. r a = yi^, and r. 
= LJLo with L„ as a representative scale length (typically a multiple of the absolute average 
deviation). The choice of scale length is a balance between a small scale for sensitivity to local 
dynamics (typically at L, < 5a) and avoidance of excessive noise (typically at L, i a). The 
distances (1^) are defined by: 

L 0 = Kr*jJ Eq. 17 

where m is the average number of points per cycle from Eq. 12 (i.e.,m = T c ). Schouten et aL (40) 
describe the details for evaluating Eqs 16- 17 to measure of the number of degrees of freedom in a 
system (e.g., the number of coupled first-order differential equations to depict the dynamics). 

The Kolmogorov entropy (K-entropy or simply entropy) measures the rate of information 
loss per unit time, or (alternatively) the degree of predictability. A positive, finite entropy 
generally is considered to be a clear demonstration that the time series and its underlying 
dynamics are chaotic. An infinite entropy indicates a stochastic, non-deterniinistic (totally 
unpredictable) phenomenon. For entropy determination, one begins with two orbits on a chaotic 
attractor that are initially very close together. The entropy then is estimated from the average 
divergence time for pairs of initiaDy-close orbits. More precisely, the entropy is obtained from the 
average time for two points on an attractor to go from an initial separation (L < L,), to become 
separated by more than a specific distance (L * L,). The maximum-likelihood entropy (K) is: 
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Eq. 19 



with b as the number of timers fit two poinu. initUny »* L<U » diverge to U W Tie 

by saouto. « * < w - — i— "«* ° f " e ^ Note *■ *" 

The entropy (K) and correlation dimension (D) usually are repotted m the limit of zero 
^ UngUt However, EEG data (and all biomedical <U») hrve subsraMnu noise. Conse,ue«ly, 
^nonhneatn^K^D^eroportedforatehesc^leng^W^is^^ 
^nrhonoise. Thus, the values of K and D, that are reported here, do not capture the full 

< «^.n■^«^|^^-- fc •-"^ fcfc, ■ M ■* 4,, •■ 

EXAMPLE I 

The tnahod »d appears of to mvamon are Dusked h, rhe of two <U»s«s. 
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designated Example L the patient was atong up m bed, domg som 
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EEG shows an activated pattern. She then showed automatisms with picking movements and 
staring, followed by vocalizations (several seconds of scream*). Hyper-extended head and neck 
posturing followed. Her upper extremities became flexed, and then she showed clonic activity, 
involving abduction/adduction at the shoulders and hips. There was tonic posturing and clonic 
activity of all extremities. The convulsive movements were associated with high-amplhude EEG 
waves, involving spikes, polyspikes, and much artifact activity. As the clinical seizur e 
spontaneously terminated, the subject was unresponsive and made loud snoring sounds. Then, the 
brain wave amplitudes became quite suppressed. The automatisms were associated with 
polyspike discharges from the left frontal region. After seizure termination, spike discharges 
occurred from this same region, followed by suppressed background waves. 

The results from Example I involved four analyses (e-, f-, g-, and h-data) on the three 
epochs of channel 13 data (epileptic seizure and post-seizure, non-seizure, and transition from 
non-seizure to seizure). These three epochs of data were provided and analyzed as three separate, 
non-contiguous ten-minute data segments. The non-seizure data segment ended several hours 
before the transition data, which in turn ended -10 seconds before the start of the seizur e data 
segment. The results are combined as one thirty-minute set of plots for the e-data (Fig. 4. 1 ), f- 
data (Fig. 4.2), g-data (Fig. 4.3), and h-data (Fig. 4.4), with gaps in the analysis to indicate where 
one dataset ends and the next dataset begins. In particular, the data from 10-590s is the non- 
seizure epoch with large amplitude signals at 520-600s from chewing and drinking. The data 
from 610-1 190s is the (transition) period immediately before seizure The data from 1210- 1790s 
includes the seizure and post-ictal phases. The various measures were obtained for every 
analysis-window, and the resulting values were plotted at the center of the 20- second analysis- 
window. 

Each plot in Figs. 4. 1 - 4.4 for Example I displays the analysis- window-centered measure 
as a solid line. The dashed line (- - -) in each figure is the average value of the measure (from Eq. 
6) over an i 1-point averaging-window, plotted at the central (sixth) point of the averaging 
window. The dotted line (...) in each figure is the corresponding sample standard deviation over 



19 



WO 97/34525 



PCT/US97/04475 



15 



this 1 l-point averagmg-window (from Eq. 9), also plotted at the central point of the averagmg- 
window. The scale length (L,) was fixed at -1.4 times the absolute average deviation, as obtamed 
by averaging over the complete non-seizure e-data, and was used in til the analyses as the 

reference scale length for all three data epochs. Smaller values for this scale length caused 

numerical problems in the determinations of the correlation dimension and the K-entropy; larger 

values limited the resolution of the nonlinear measures. 

The clinical seizure in Example I occurred from 53 to 95s in the (third) seizure epoch 

(1253 to 1295sinFigs.4.M.4). Rhythmic convulsions began at 1295s, and post-ictal features 

appeared at 1314s. 

T„nds m the average and deviate. (SD) of each measure forth* dam*.™ 

very diHerou. during <»>«-l »*> " — P— - *• ^ 
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crroWon dimeutno. . C) The SD of each tneuuro (x) is donored by (e.g., -he SD of.be 
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trend also occurs in K and Q of g-data at 690s (maxX 790s (min), 880s (max), 965s (min), and 
1060s (max), with corresponding minima and maxima in £ for the first four times. This trend 
occurs again in & and D of h-data at 695s (max), 790s (min), 870s (max), 960s (min), and 1060s 
(max), with corresponding minima and maxima in £ that approximate the first four times. The 
times at which the cxtrema occur is within ±10s among to e-/g-/h-data, and is consistent with 
uncertainty that is inherent in the 10s overlap of successive analysis windows. 

A third trend in the g-data of Example I involves a valley in M t that is coincident with a 
peak in o M . The rise in o M begins at 785s, peaking at 885s, with the subsequent decrease ending at 
990s. The corresponding times for Mi (i.e., beginning of the falL minimum, and end of the rise) 
are coincident with the corresponding time values for o M . This trend in M, and o M corresponds to 
oscillations in M, that increase in amplitude from 785-885s and subsequently decay from 885- 
990s. Note that the extrema in M, and o,, occur at 885s, which is approximately coincident with 
one of the extrema times in (K . D. TJ as discussed in the previous paragraph. 

Other trends are inferred in Example I from the total variation in the average and in the 
standard deviation of each measure during the non-seizure period relative to that during the 
transition period. For example, the average skewedness in e-data (see Fig. 4. 1 ) is -0. 52 i 5 s 
0.25 (during the non-seizure period) and -0.5 1 s s s 0.24 (during transition). The resulting total 
variations (A) in average are A, = 0.77 (with the subscript "n" denoting non-seizure) and A, = 
0.75 (with the subscript "t" denoting transition). The ratio of these variations is A/A„ = 0.97. 
Likewise, the total variation in the SD of the skewedness of e-data is Ao B =0.9 (non-seizure) and 
Ao, = 0.6 (transition). Then, the ratio of the total variation in standard deviation is Ao/Ao, = 
0.67. Table 1 summarizes these comparisons between non-seizure and transition data, taken from 
Figs. 4. 1-4.4. Most of these variations are not significant because noise in the EEG data (and 
biological data in general) makes comparisons tenuous for ratios near unity. However, starred 
values (*) in Table 1 are probably significant, corresponding to substantial ratios (e.g. >2 or 
<0.5). We note that some small ratios (sO.l) exist only for the linear measures of g-data (e.g., 
A/A 0 = 0.016 and Ao/Ao B = 0.082 for kurtosis, and Ao/Ao a = 0. 1 1 for skewedness). In this 
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patient, these variations mean that the non-seizure EEG is much more variable than the transition 
EEG for g-data for this type of seizure, and the transition is detectable with linear and non-linear 
statistics. 

Listed below are indicators of transition from non-seizure to seizure found in Example I: 
O, i 0.5 for J-200 sec. for e-data and f-data 
quasi-periodic mhVmax in K of e-data, g-data, and h-data 
quasi-periodic mm/max in E* of e-data, g-data, and h-data 
quasi-periodic mm/max in I.of e-data, g-data, and h-data 
valley in Mi coincident with peak in o m of g-data 
A/A, s 0.5 for g-data (skewedness) 
Ao/Ao B < 0.5 for f-datt and g-data (skewedness) 
A/A, s 0.5 for e-data, g-data, and h-data (kurtosis) 
Ao/Ao n < 0.5 for e-data, g-data, and h-data (kurtosis) 
A/A„ < 0.5 for g-data (T.) 
Ao,/Ao n <0.5 for g-data and h-data (TJ 
Ao/Ao^O.5 for f-data and h-data (K) 
A/A n <0.5 for f-data (M,) 
Ao/Ao^O.5 for f-data (M,) 
A/A,<0.5 for g-data (D) 
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Table 1 Ratios of averages and standard deviations for Example I 



Specific measure e-data f-data g-data h-dau 



skewedness (s) 

A/A, 0.97 1.1 0.29* 0.77 

Ao/Ao B 0.67 0.48* 0.11* 0.73 

kurtosis (k) ^ 

A/A. 0.37* 0.90 0.016* 0.44* 

Ao/Ao n 0.28* 0.58 0.082* 0.42* 



Tone per cycle (T c ) 

A/A. 0.51 0.67 0.35* 0.63 



Ao/Ao B 0.67 0.62 0.35* 0.46 
Entropy (K) 

A/A„ 0.71 0.50 0.70 0.53 

Ao/Ao, 0.93 0.33* 0.82 0.42* 

rMin. inMIFCM,) 

A/A, 1.3 0.40* 0.52 0.67 

Ao/Ao, 0.51 0.40* 0.91 0.50 

Correlation dimension (D) 

A/A„ 0.56 0.73 0.39* 0.97 

Ao/Ao, 0.64 0.57 0.50 0.78 
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EXAMPLE 0 

Hie dataset designated Example D is from the same patient as Example L who is a 20- 
year-old female with a lifelong history of seizures as previously described m Example I. During 
the seizure designated Example n, the patient is tying in bed, awake with her right upper 
extremity in a flexed posture. The EEG shows spike-wave discharges in the left hemisphere, after 
which the brain waves become very sharp, dominated by high frequency activity and artifacts. 
The patient showed eye deviation to the right and some head turning to the right, with head 
jerking also to the right, but without any usual posturing of the right upper extremity. Thc 
eye/head turning is preceded by a high-frequency vocalization. After the clinical seizure 
spontaneously terminated, the EEG shows high amplitude wave slowing and subsequent 
amplitude suppression. The subject remained awake during the seizure but was poorly 
responsive. 

The results from Example H for channel 13 involved four analyses (e-, f-, g-, and h-dau) 
on one 23-mmute dataset that included all three epochs of channel 13 data (epileptic seizure, non- 
seizure andtransmonfromnon-seiziiretoseizure). The non-seizure period spanned IMOOs. 
The transition period occurred over 410-1200*. The seizure began at 1245s and ended at 1290s. . 
The patient was aphasic at 1300s, with head movements and verbalization at 1315s. 

The scale length (L,) was fixed at -1.0 times the absolute average deviation, as obtained 
by averaging over the non-seizure e-data, and was used in all the analyses as the reference scale 
length for all three data epochs. The results were obtained as before, and are plotted for e-data 
(Rg. 5.1), f-data, (Fig. 5.2), g-data (Rg. 5.3), andh-data (fig. 5.4). 

As in Example I, certain trends in the data are useful as seizure predictors for this pauent 
One such trend in Example U involves the sample standard deviation of the correlation dimension 
(o D ) for e-h data. In particular, o D has a small variation during the non-seizure period, but 
undergoes large rises and fans during transition. Moreover, o D has minimal variation (o D , 0.1) 
for 30-150S intervals during transition. This trend is unique to the transition period. 

A second trend in Example H involves the entropy (K) for f- and h-data. In these cases, K 
has a large, aperiodic variations during transition, and substantially less variation during the non- 
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seizure period Also during transition, the entropy becomes very small (K < 0.00005) for several 
intervals that have a duration of 10- 100s. 

A third trend in Example II occurs in the first minimum of the MIF (Mj). From 700- 
1100s, Mt in the e-data gradually decreases below a nominal critical value (40 in this example) to 

S a minimum of 21 just before the seizure. A similar decrease occurs in Mi from 400-1 100s in g- 
data, to a mnrimiim of 17. Two cycles of a quasi-periodic variation of Mi in f-data appear from 
4OO-800s, with a valley-to-valley period of 200s. 

A fourth trend in Example II involves the entropy (K) for e-data. From 600-730s, the 
average entropy (K) increases monotonically, while K undergoes oscillations of decreasing 

10 amplitude about K. The entropy then decreases abruptly (K=0.045 -> K==0.0004) during 760- 
820s. Thai, K again increases monotonically over 85 0-920 s, while K undergoes oscillations of 
increasing amplitude about this average. The entropy then decreases abruptly (K=0.0 18 — 
10=0.001) during 930-970s. The previous cycle is repeated over 1000- 1 100s (increasing K with K 
undergoing oscillations of increasing amplitude), followed by a smaller decrease in (K=0.009 

IS K=0.003). This last decrease is difficult to distinguish from the imbedded oscillations, which 
subsequently rise again as part of the seizure. 

Other trends in the average and standard deviation of each measure for Example II are 
very different during the transition period (4 10- 1200s) as compared to the non-seizure period ( 10- 
400s). Table 2 shows the ratios of the variations in the average and in the standard deviation of 

20 each measure during the non-seizure period with that during the transition period. As before, 
most of these variations are not significant because noise in the EEG data makes comparisons 
tenuous for ratios near unity. However, starred values (*) in Table 2 are probably significant, 
corresponding to substantial ratios (e.g. >2 or <0.5). Note that some large ratios (>10) occur for 
the linear measures in Table 2. These large ratios mean that the transition data is much more 

25 variable for this seizure than for the non-seizure data. 

Listed below are indicators of transition from non-seizure to seizure found in Example II: 
(o D ) D has small variation, (o D X has large variation (e-data, f-data, g-data, and h-data) 
(K\ has large variation,^), has much less variation (f-data and h-data) 
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(KX becomes very small (< 5x10 5 ) for 10 • 100s 
gradual decrease in Mi (e-data and g-data) 
quasi-periodic variation in M, (f-data) 

K increases monotonically while K oscillates about £ (e-data), followed by an abrupt decrease in 
5 K 

A/A, > 2 for SD of e-data, f-data, g-data, and h-data 
Ao/Ao„ > 2 for SD of e-data, f-data, and h-data 
A/A, > 2 for AAD of e-data, f-data, g-data, and h-data 
Ao/Ao, > 2 for AAD of e-data, f-data, and h-data 
10 A/A,, > 2 for skewedness of e-dau, f-data, g-data, and h-data 

Ao/Ao n > 2 for skewedness of e-data, f-data, g-data, and h-data 

A/A, > 2 for kurtosis of e-data, f-data, g-data, and h-data 

Ao/Ao, > 2 for kurtosis of e-dau, f-data, g-data, and h-data 

Ao/Ao,< 0.5 for K of g-data 
15 A/A. >- 2 for M t off-data, g-data, and h-data 

A/A, > 2 for D of h-data 

A/A, v 2 for D of e-data and h-data. 
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Table 2 Ratios of averages and standard deviations for Example U 
Specific measure e-data f.data g-data h-data 

std. deviation (o) 

A/A, 3.9* 4.8* 2.1* 3.3* 

Ao/Ao B 15.7* 25.9* 1.2 6.5* 



abs. avg. dev. (a) 

A/A„ 2.5* 2.9* 2.3* 3.0* 

Ao/Ao. 10.6* 12.7* i.i 3 . 6 . 

skewedness (s) 

A/A. 4.8* 4.1* 8.6* 5.5* 



Ao/Ao, 3.1* 3.8* 



44.5* 23.9' 



kurtosis (k) 

6.9* 8.6* 10.8* 17.5 

Ao/Ao„ ii.7* 69 . 36 ,. M2 



18 1.0 
1.0 0.67 



Time per cycle (TJ 

A/A, 1.4 1.3 

Ao/Ao, 1.2 i.g 

Entropy (K) 

^ 1-5 1.9 1.00 2.0 

Ao/Ao B 0.99 1.5 0.49* 1.6 
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Table 2 Continued 

l»Min.inMIF(M,) 

A/A„ " 30 * 26 * 35 * 



1.4 1.1 1-7 2° 



5 Correlation dimension (D) 
A/A. 
Ao/Ao B 



1.8 
4.3* 



1.8 
1.9 



1.7 
1.7 



3.3* 
2.2* 



10 



15 



20 



25 



It has been shown by the foregoing that certain trends such as sudden increases, sudden 
decreases, peaks, and valleys and combinations thereof in time serial sequences of nonlinear 
measures hereinbefore described are seizure predictors. 

By comparing indicative trends in non-linear measures of the patient's brain wave data 
with known seizure indicators, it can be clearly determined whether the indicative trends 
correspond with known seizure predictors whether a seizure is oncoming in the patient. 

Tables 1 and 2 and Figures 4 and 5 clearly indicate that epileptic seizures are mdicated by 
a number of trends in nonlinear measures of brain-wave data mcluding the following: 

a a long-time, large value of the standard deviation of the skewedness of e-data 101; 
b a long-time, large value of the standard deviation of the skewedness off-data 103; 
c quasi-periodic maxima and minima in the average entropy of e-data 105; 
d quasi-periodic maxima and minima in the average entropy of g-data 107; 

e quasi-periodic maxima and minima in the average entropy of h-data 109; 

f quasi-periodic maxima and minima in the average correlation dimension of e-data 111 . 

g quasi-periodic maxima and nmumain the avenge conebtion dimension ofg-datall3; 

h quasi-periodic maxima and minima in the average correlation dimension of h-data 115; 

i quasi-periodic maxima and minima in the average 1 
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j quasi-periodic maxima and minima in the average number of time steps per cycle of g-data 119; 
k quasi-periodic maxima and minima in the average number of time steps per cycle of h-data 121; 
1 coincidence of the minima and maxima selected from the group consisting of items c) 

through k) of this list and combinations thereof; 
m a valley in the average of the minimum m the mutual information function in g-data 122; 
n a peak m the standard deviation of the minimum m the mutual information function in g-data 

123; 

o coincidence of the extrema selected from the group consisting of hems m) and n) of this list 
and combinations thereof with at least one of extrema times selected from the group 
consisting of hems c) through k) of this list and combinations thereof; 

p large decreases and increases in the standard deviation of the correlation dimension of e-data 
l25 > 

q large decreases and increases in the standard deviation of the correlation dimension off-data 
127; 

r large decreases and increases in the standard deviation of the correlation dimension of g-data 
129; 

s large decreases and increases in the standard deviation of the correlation dimension of h-data 
131; 

t large, aperiodic variations in entropy of f-data, including some values close to zero for > 10 
seconds 133; 

u large, aperiodic variations in entropy of h-data, including some values close to zero for >10 
seconds 135; 

v the average of the first minimum of the mutual information function of e-data gradually 

decreases below a critical vahiel37; 
w the average of the first nunimum of the mutual information function of g-data gradually 

decreases below a critical value 139; 
x quasi-periodic maxima and minima in the average of the first imnimum of the mutual 

information function off-data 141; 
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y aperiodic maxima and minima in the average entropy of e-data 143; 

z escalations in entropy of e-data of increasing or decreasing magnitude, about the average 

entropy of item y 145; 
aa small variation ratio in average skewedness (A/A.) of g-data 147; 
ab small variation ratio in standard deviation (Ao/AoJ of skewedness of f-datal49; 
ac small variation ratio in standard deviation (Ao/Ao.) of skewedness of g-data 151; 
ad small variation ratio in average kurtosis (A/A.) of e-data 153; 
ae small variation ratio in average kurtosis (A/A„) of g-data 155; 
af small variation ratio in average kurtosis (A/AJ of h-data 157; 
ag small variation ratio in standard deviation (Ao/AoJ of kurtosis of e-data 159; 
ah small variation ratio in standard deviation (Ao/AoJ of kurtosis of g-data 161; 
ai small variation ratio in standard deviation (Ao/AoJ of kurtosis of h-data 163; 
aj smanvariationnmomavengenu^^ 
• ak small variation ratio in standard deviation (Ao/Ao J of number of timesteps per cycle of g- 

data 167; 

a small variation ratio in standard deviation (Ao/Ao,) of nornberoftinKsaeps per cycle of b- 
data 169; 

am small variation ratio in standard deviation (Ao/Ao,) of entropy of f-datal71; 
aa smaEvariarionrariomstandarddeviationCAo./AoJof entropyof h-data 173; 
ao smanvariauonratmmmeaverageofme^m^ 
of f-data 175; 

ap small variation ratio in standard deviation (Ao/AoJ of Srssmrainmmmnmro.lmforma.ron 

function off-data 177; 
aq small variation ratio in average correlation dimension (A/A n ) of g-data 179; 
ar large variation ratio in average (A/A.) of the standard deviation of e-data 181 ; 

~ «,*roo* /A7A ^ of the standard deviation of f-data 183; 
as large variation ratio m average 01 nauuui 

. . „ fimvP (\JK\ of the standard deviation of g-data 185; 
at large variation ratio m average 01 ^ 
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au large variation ratio in average (A/AJ of the standard deviation of h-data 187; 
av large variation ratio in standard deviation (Ao/Ao J of the standard deviation of e-data 189; 
aw large variation ratio in standard deviation (Ao/AoJ of the standard deviation of f-data 191 
ax large variation ratio in standard deviation (Ao/AdJ of the standard deviation of h-data 193: 
ay large variation ratio in the average of the absolute average deviation (A/A,) of e-data 195; 
az large variation ratio m the average of the absolute average deviation (A/AJ of f-data 197; 
ba large variation ratio in the average of the absolute average deviation (A/A,) of g-data 199; 
bb large variation ratio in the average of the absolute average deviation (A/AJ of h-data 201; 
be large variation ratio in the standard deviation of the absolute average deviation (Ao/AoJ of 
e-data 203; 

bd large variation ratio in the standard deviation of the absolute average deviation (Ao/AoJ of 
f-data 205; 

be large variation ratio in the standard deviation of the absolute average deviation (Ao/Ao J of 
h-data 207; 

bf large variation ratio in average skewedness (A/AJ of e-data 209; 

bg large variation ratio in average skewedness (A/AJ of f-data 211; 

bh large variation ratio in average skewedness (A/A.) of g-data 213; 

bi large variation ratio in average skewedness (A/AJ of h-data 215; 

bj large variation ratio in standard deviation (Ao/AoJ of skewedness of e-data 217; 

bk large variation ratio in standard deviation (Ao/Ao J of skewedness of f-data 219; 

bl large variation ratio in standard deviation (Ao/Ao,) of skewedness of g-data 221; 

bm large variation ratio in standard deviation (Ao/Ao J of skewedness of h-data 223; 

on large variation ratio m average kurtosis (A/AJ cf e-data 225; 

bo large variation ratio m average kurtosis (A/A.) of f-data 227; 

bp large variation ratio in average kurtosis (A/AJ of g-data 229; 

bq large variation ratio in average kurtosis (A/AJ of h-data 231 ; 

br large variation ratio in standard deviation (A0/A0J of kurtosis of e-data 233; 
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bs large variation ratio m standard deviation (Ao/Ao.) of kurtosis of f-data 235; 
be large variation ratio in standard deviation (Ao/Ao n ) of kurtosis of g-data 237; 
bu large variation ratio in standard deviation (Ao/Ao B ) of kurtosis of h-data 239; 
bv small variation ratio in standard deviation (Aa/Ao,) of entropy of g-data 241; 
5 bw large variation ratio in the average of the first nriniznum in mutual information function 
(A/AJ of f-data 243; 

bx large variation ratio in the average of the first nunimum in mutual mformauon function (A/A.) 

by Lr^r^ 

10 of h-data 247; 

bz tote v«i^ ratio « cond.do. (A/W of W«. 2«; 

« ^varato ios^ deviation (Ao/Ao.) of a^ral^ dimannon of 

* 1 ^ e v M Woox«io to su.d M ddov i a t io.(Ao/AoJof conoUuon dtaaaion of Man. 253; 

and 

15 cc combinations thereof 

White thara hasbaen shown Bid described what - pre»Kly coosdarad to prdanad 
- ^^l-^*-*.«W--*-**» , '*'"**" ,, ".*71 1 

by the appended claims. 
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What is claimed is: 



(1) 



A method for automatically predicting an epileptic seizure in a patient comprising the steps 
of 

(a) providing at least one channel of a patient's raw brain wave data, called e-data, 
selected from the group consisting of electroencephalogram data and 
magnetoencephalogram data; 

(b) separating the e-data into artifact data, called f-data, and artifact-free data, called 
g-data, while preventing phase distortions in the data; 

(c) processing g-data through a low-pass filter to produce a low-pass-filtered version 
of g-data, called h-data; 

(d) applying at least one measure selected from the group consisting of the linear 
statistical measures mhumum and maximum, standard deviation, absolute minimum, 
deviation, skewedness, and kurtosis, and the nonlinear measures time steps per 
cycle, Kolmogorov entropy, first minimum in the mutual information function, and 
correlation dimension to at least one type of data selected from the group 
consisting of e-data, f-data, g-data, and h-data to provide at least one time serial 
sequence of nonlinear measures, from which at least one indicative trend selected 
from the group consisting of abrupt increases, abrupt decreases, peaks, valleys, 
and combinations thereof is determined; 

comparing at least one indicative trend with at least one known seizure predictor; 
and 



(c) 
(f) 



determining from said comparison whether an epileptic seizure is oncoming in the 
patient 
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(2) 



The method as described in Claim 1 wherein said at least one time serial sequence oflinear 
and nonlinear nonlinear measures is selected from the group consisting of. the sample 
standard deviation of the correlation dimension for e-data; the sample standard deviation 
of the conelation dimension for f-data; the sample standard deviation of the conelation 
dimension for g-data; the sample standard deviation of the conelation dimension for h- 
data; the Kolmogorov entropy for f-data; the Kohnogorov entropy for h-data; the first 
mutual to the Mutual Information Function for e-data; the first mutual to the Mutual 
Information Function for g-data; the average Kolmogorov entropy for e-data; the 
skewedness for f-data; the kurtosis for f-data, and combinations thereof 

(3) The method as described to Claim 2 whereto said at least one seizure predictor is selected 
from the group consisting of. 

(a) a long-time, large value of the standard deviation of the skewedness of e-data; 

(b) a long-time, large value of the standard deviation of the skewedness off-data; 
quasi-periodic maxima and imnima to the average enttopy of e-data; 
quasi-periodic maxima and minima to the average entropy of g-data; 
quasi-periodic maxima and minima to the average entropy of h-data; 
quasi-periodic maxima and ininima to the average conelation dimension ofc- 
data; 

quasi-periodic maxima and minima to the average conelation dimension ofg- 
data; 

quasi-periodic maxima and minima to the average correlation dimension ofh- 
data; 

quasi-periodic maxima and minima to the average number of time steps per cycle 
of e-data; 

quasi-periodic maxima and minima to the average number of time steps per cycle 
ofg-data; 

quasi-periodic maxima and imnima to the average number of time steps per cycle 



(c) 
(d) 
(c) 
(f) 



(g) 
00 
(i) 
(j) 
00 
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ofh-data; 

(1) coincidence of the minima and maxima selected from the group consisting of 

items c through k) of this list and combinations thereof; 
(m) a valley in the average of the minimum in the mutual information function in g- 

data; 

(n) a peak m the standard deviation of the minimum in the mutual information 
function in g-data; 

(o) coincidence of the extrema selected from the group consisting of items m and 

n of this list and combinations thereof with at least one of extrema times selected 
from the group consisting of items c through k of this list and combinations 
thereof; 

(p) large decreases and increases in the standard deviation of the correlation 
dimension of e- data; 

(q) large decreases and increases in the standard deviation of the correlation 
dimension of f-data; 

(r) large decreases and increases in the standard deviation of the correlation 
dimension of g-data; 

(s) large decreases and increases in the standard deviation of the correlation 
dimension ofh-data; 

(t) large, aperiodic variations in entropy of f-data, including some values close to 

zero for >10 seconds; 
(u) large, aperiodic variations in entropy of h-data, including some values close to 

zero for > 10 seconds; 
(v) the average of the first minimum of the mutual information function of e-data 

gradually decreases below a critical value; 
(w) the average of the first minnnnm of the mutual information function of g-data 

gradually decreases below a critical value; 
(x) quasi-periodic maxima and minima in the average of the first mmiim^ 0 f ^ 
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mutual information function off-data; 
(y) aperiodic maxima and minima in the average entropy of e-data; 
(z) oscillations in entropy of e-data of increasing or decreasing magnitude, about the 

average entropy of item y; 

(aa) small variation ratio m average skewedness (A/AJ of g-data; 

(ab) smaflvariationntiomstandarddev*^^ 

(ac) small variation ratio in standard deviation (Ao/AoJ of skewedness of g-dau; 

(ad) small variation ratio in average kurtosis (A/AJ of e-dau; 

(ae) small variation ratio in average kurtosis (A/A.) of g-data; 

(af) small variation ratio in average kurtosis (A/A„) of h-data; 

(ag) small variation ratio in standard deviation (Ao/Ao„) of kurtosis of e-data; 

(ah) smaU variation ratio in standard deviation (Ao/Ao B ) of kurtosis of g-data; 

( ai) small variation ratio in standard deviation (Ao/AoJ of kurtosis of h-data; 
(aj) small variation ratio in average number of timesteps per cycle (A/AJ of g-data; 
(ak) small variation ratio in standard deviation <Ao/Ao B ) of number of timesteps per 

cycle of g-data; 

(al) . small variation ratio in standard deviation (Ao/Ao n ) of number of timesteps per 
cycle of b-data; 

(am) small variation ratio in standard deviation (Ao/AoJ of entropy of f-data; 
(an) small variation ratio in standard deviation (Ao/Ao.) of entropy of h-data; 
(ao) small vanation ratio in the average of the first minimum m mutual informal 

fimction(A/AJ of f-data; 
(ap) small variation ratio in standard deviation ( Ao/Ao.) of first mhtimum m mutual 

information function of f-data; 
(aq) small variation ratio in average correlation dimension <A/A n ) of g-dau; 
(ar) large variation ratio in average (A/A,) of the standard deviation of e-data; 
(as) large variation ratio in average (A/A,) of the standard deviation of f-data; 
(at) large variation ratio in average (A/AJ of the standard deviation of g-dau; 
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(au) large variation ratio in average (A/AJ of the standard deviation of h-data; 
(av) large variation ratio in standard deviation (Ao/Ao.) of the standard deviation of 
e-data; 

(aw) large variation ratio in standard deviation (Ao/Ao.) of the standard deviation of 
f-data; 

(ax) large variation ratio in standard deviation (Ao/AoJ of the standard deviation of 
h-data; 

(ay) large variation ratio in the average of the absolute average deviation (A/A.) of e- 
data; 

(az) large variation ratio in the average of the absolute average deviation (A/AJ of f- 
data; 

(ba) large variation ratio in the average of the absolute average deviation (A/AJ of g- 
data; 

(bb) large variation ratio in the average of the absolute average deviation (A/A.) of fa- 
data; 

(be) large variation ratio in the standard deviation of the absolute average deviation 
(Ao/Ao.) of e-data; 

(bd) large variation ratio in the standard deviation of the absolute average deviation 
(Ao/Ao.) off-data; 

(be) large variation ratio in the standard deviation of the absolute average deviation 
(A0/A0J of h-data; 

(bf) large variation ratio in average skewedness (A/A.) of e-data; 

(bg) large variation ratio in average skewedness (A/A.) of f-data; 

(bh) large variation ratio in average skewedness (A/A.) of g-data; 

(bi) large variation ratio in average skewedness (A/AJ of h-data; 

(bj) large variation ratio in standard deviation (Ao/Ao.) of skewedness of e-data; 
(bk) large variation ratio in standard deviation (Ao/Ao.) of skewedness of f-data; 
(bl) large variation ratio in standard deviation (Ao/Ao.) of skewedness of g-data; 
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(bin) large variation ratio in standard deviation (Ao/Ao.) of skewedness of h-data; 

(bn) large variation ratio in average kurtosis (A/A B ) of e-data; 

105 (bo) large variation ratio in average kurtosis (A/A B ) of f-data; 

(bp) large variation ratio in average kurtosis (A/A,) of g-data; 

(bq) large variation ratio in average kurtosis (A/A.) of h-data; 

(b r) large variation ratio in standard deviation (Ao/Ao B ) of kurtosis of e-data; 

(bs) large variation ratio in standard deviation (Ao/Ao.) of kurtosis of f-data; 

(bt) large variation ratio in standard deviation (Ao/Ao.) of kurtosis of g-data; 

(bu) ^v^r**<^*^WW«1^«^ 

(bv) sn uUvariauonratiomstandarddcviation(Ao/Ao. ) of entropy of g-data; 

(bw) l BI ,^ii4.«^rf* ta ^ k,rt,d,,-k, 
&nction(A,/AJ of f-data; 

• • .u- .„ mo * of the first minimum in mutual information 
(b x ) large variation ratio m the average ottnenrsi minimum 

function (A/AJ of g-data; 

• «„n ™ti« in the average of the first mhumnm in mutual infonnauoo 
(by) large vanatton ratio m tne average ui u»*. 

function (A/A«) of h-data; 



(bz) large variation ratio in average correlation dimension (A/AJ of h-data^ 
e-data; 



large vjuwuuu ~ o 

e-data; f 



h-data; and 
(cc) combinations thereof 

(4) U e ^od K d^ to C U ml^^«i & »^^ mte4& ° m,ht,,W 
data by use of a zero-phase finer. 

, . - . w, Claim 1 wherein said low-pass filter comprises a standard low- 

(5) The method as described m Claim l wnnnu ^ r 

y . ^«» of second-order, third-order and fourth- 

pass filter selected from the group consisting of second outer, uuru 
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order low-pass filters at frequencies between about 35 Hz and about 60 hi 

The method as described in Claim 5 wherein said low-pass filter comprises a standard 
fourth-order low-pass filter at 50 Hz. 

Apparatus for automatically predicting an epileptic seizure in a patient comprising: 

(a) data provision means for providing at least one channel of raw brain wave data, 
caUed e-data, selected from the group consisting of electroencephalogram data and 
magneto encephalogram data: 

(b) separation means for separating e-data into artifact data, called f-data, and artifact- 
free data, called g-data. while preventing phase distortions in the data, 
communicably connected to said data provision m eans ; 

(c) low-pass filter means for filtering g-data to produce a low-pass filtered version of g- 
data, called h-data, communicably connected to said separation means: 

(d) application means for applying at least one measure selected from the group of 
consisting of the linear statistical measures minimum and maximum, standard 
deviation, absolute minimum deviation, skewedness, and kuitosis, and the 
nonlinear measures time steps per cycle, kolmogorov entropy, first minimiTm in the 
mutual information function, and correlation dimension to at least one type of data 
selected from the group consisting of e-data, f-data, g-data, and h-data to provide at 
least one time serial sequence of nonlinear measures, from which at least one 
indicative trend selected from the group consisting of abrupt increases, abrupt 
decreases, peaks, valleys, and combinations thereof is determined, communicably 
connected to said low-pass filter means; 

(e) comparison means for comparing at least one indicative trend with known seizure 
predictors, connected to said application means; and, 

(f) determination means for determining from the comparison whether an epileptic 
seizure is oncoming in the patient communicably connected to said comparison 
means. 
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(8) The apparatus as described in Claim 7 wherein said at least one time serial sequence of 
linear and nonlinear nonlinear measures is selerted from me group consisting of: the 
sample standard deviation of the correlation dimension for e-data; the sample standard 
deviation of the correlation dimension for f-data; the sample standard deviation of the 
correlation dimension for g-data; the sample standard deviation of the correlation 
dimension for h-data; the Kolmogorov entropy for f-data; the Kolmogorov entropy for h- 
data; the first mutual in the Mutual Information Function for e-data; the first mutual in the 
Mutual Information Function for g-data; the average Kolmogorov entropy for e-data; the 
skewedness for f-data; the kurtosis for f-data, and combinations thereof 

(9) The apparatus as described in Claim 8 wherein said at least one seizure predictor is 
selected from the group consisting of 

a long-time, large value of the standard deviation of the skewedness of e-data; 
a long-time, large value of the standard deviation of the skewedness off-data; 
quasi-periodic maxima and ininima in the average entropy of e-data; 
quasi-periodic maxima and minima in the average entropy of g-data; 
quasi-periodic maxima and minima in the average entropy of h-data; 
quasi-periodic maxima and minima m me average conebtion dimension ofe- 

data; 

quasi-periodic maxima and minima in the average correlation dimension of g- 
data; 

quasi-pcriodic maxima and minima in the average correlation dimension of h- 
data; 

quasi-periodic maxima and minima in the average number of time steps per cycle 
of e-data; 

j m the averaee number of time steps per cycle 

quasi-periodic maxima and minima m tne average uumu 

ofg-data; 



(») 
(b) 

(c) 
(d) 
(c) 
(f) 

(g) 

(h) 

(0 

0) 
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(k) quasi-periodic maxima and minima in the average number of time steps per cycle 
ofh-data; 

(1) coincidence of the minima and maxima selected from the group consisting of 

items c through k of this list and combinations thereof; 
(m) a valley in the average of the minimum in the mutual information function in g- 

data; 

(n) a peak in the standard deviation of the minimum in the mutual information 
function in g-data; 

(o) coincidence of the extrema selected from the group consisting of hems m and 
n of this list and combinations thereof with at least one of extrema times selected 
from the group consisting of hems c through k of this fist and combinations 
thereof; 

(p) large decreases and increases in the standard deviation of the correlation 
dimension of e-data; 

(q) large decreases and increases in the standard deviation of the correlation 
dimension off-data; 

(r) large decreases and increases in the standard deviation of the correlation 
dimension of g-data; 

(s) large decreases and increases in the standard deviation of the correlation 
dimension ofh-data; 

(t) large, aperiodic variations in entropy of f-data, including some values close to 

zero for > 10 seconds; 
(u) large, aperiodic variations in entropy of h-data, including some values close to 

zero for >10 seconds; 
(v) the average of the first nrinimum of the mutual information function of e-data 

gradually decreases below a critical value; 
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(w) 

gradually decreases below a critical value; 

a — u.;^. ' m the average of the first mnritTnitn of the 
(x) quasi-periodic maxima and mmnnam the avenge m« 

mutual information function of f-data; 
(v) aperiodic maxima and nimiinam the average entropy of e-data; 
Z oscillations in entropy of e-data of tncreasing or decreeing magnitude, about the 

average entropy of item y; 
(aa) small variation ratio in average skewedness (A/AJ of g-data; 

(K ) ^ variation ratio in standard deviation (Ao/AoJ of skewedness of g-data; 
(ad) small variation ratio in average kurtosis (A/AJ of **»*; 
( ac) smanvariauonratiomaveragekurtosisCA/^ofg-data; 



(af) small variation ratio in average knrtosis (A,/ A.) of h-data; 

) ■ 
) 
) 
) 

0 

cycle of g-data; 



J 

. j I AoJAo > of number of nmcaeps pa 

(t k) small vaiiition ratio m sumd«rd devuraoB (floyoo.) 01 



(ai) smtf v^rauom^dev^Ao/AoJof tooasof 

' ^ va ^o.ranomsraod^d=vuto(Ao^Jof tan .sis.f^; 
^v^nooratiora^raoOoun^of^epspercycMA/AJof^ra; 



(al) Inv.n^.ran.rasrandudd^tlooWAoJof o»mb e rof 1 ime»cpsp» 
cycle of h-data; 

(a m) aMB «n.^ra ti orasra,d M dd W rauo 1 (Ao^.f ^ • ««• 
<J> s^van^ra^ra.h.^.ofraotoo^raraunra.mf.n^ 
taata.(A/AJof f-d«u; , c .. ..^j 

information function off-data; 
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(aq) small variation ratio in average correlation dimension (A/A,) of g-data; 
(ar) large variation ratio in average (A/AJ of the standard deviation of e-data; 
(as) large variation ratio in average (A/AJ of the standard deviation of f-data; 
(at) large variation ratio in average (A/AJ of the standard deviation of g-data; 
(au) large variation ratio in average (A/AJ of the standard deviation of h-data; 
(av) large variation ratio in standard deviation (AO/A0J of the standard deviation of 
e-data; 

(aw) large variation ratio in standard deviation (Ao/AoJ of the standard deviation of 
f-data; 

(ax) large variation ratio in standard deviation (A0/A0J of the standard deviation of 
h-data; 

(ay) large variation ratio in the average of the absolute average deviation (A/AJ of e- 
data; 

(az) large variation ratio in the average of the absolute average deviation (A/AJ of f- 
data; 

(ba) large variation ratio in the average of the absolute average deviation (A/AJ of g- 
data; 

(bb) large variation ratio in the average of the absolute average deviation (A/AJ of fa- 
data; 

(be) large variation ratio in the standard deviation of the absolute average deviation 
(A0/A0J of e-data; 

(bd) large variation ratio in the standard deviation of the absolute average deviation 
(A0/A0J off-data; 

(be) large variation ratio in the standard deviation of the absolute average deviation 
(A0/A0J of h-data; 

(bf) large variation ratio in average skewedness (A/A,) of e-data; 

(bg) large variation ratio in average skewedness (A/AJ of f-data; 

(bh) large variation ratio in average skewedness (A/AJ of g-data; 
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(bi) large variation ratio in average skewedness (A/A B ) ot h-data; 

(bj) large variation ratio in standard deviation (Ao/Ao.) of skewedness of e-data; 

(bk) large variation ratio in standard deviation (Ao/Ao.) of skewedness of f-data; 

(bl) large variation ratio in standard deviation (Ao/Ao 0 ) of skewedness of g-data; 

(bm) large variation ratio in standard deviation (Ao/Ao.) of skewedness of h-data; 

(bn) large variation ratio in average kurtosis (A/A,) of e-data; 

(bo) large variation ratio in average kurtosis (A/A„) of f-data; 

(bp) large variation ratio in average kurtosis (A/AJ of g-data; 

(bq) large variation ratio in average kurtosis (A/A„) of h-data; 

(br) large variation ratio in standard deviation (Ao/Ao.) of kurtosis of e-dau; 

(bs) large variation ratio in standard deviation (A0/A0J of kurtosis of f-data; 

(bt) large variation ratio in standard deviation (Ao/Ao.) of kurtosis of g-data; 

(bu) large variation ratio in standard deviation (A0/A0J of kurtosis of h-data; 

(bv) small variation ratio in standard deviation (Ao/Ao.) of entropy of g-data; 

(bw) large variation ratio in the average of the first irnn^ 

function (A/A.) of f-data; 
(bx) large variation ratio m me average of me fu^ 

function (A/A.) of g-data; 
(by) large variation ratio in the average of the fcst nunirnum in nmtual Monnation 

function (A/AJ of h-data; 
(bz) large variation ratio in average correlation dimension (A/A.) of h-data; 
(c.) large variation ratio in standard deviation ( Ao/Ao.) of correlation dimension of 

e-data; 

(cb) large variation ratio in standard deviation (Ao/Ao.) of correlation dimension of 

h-data; and 

(cc) combinations thereof 

(10) Tne apparatus as described in Claim 7 wherein said separation means comprises a zerc- 
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phase filter. 

(11) The apparatus as described in Claim 10 wherein said zero-phase filter is embodied in a 
programmed integrated-circuit semiconductor chip. 

( 12) The apparatus as described in Claim 7 wherein said low-pass fiber means comprises a 
standard low -pass filter selected from the group consisting of second-order, third-order, 
and fourth-order low-pass filters at frequencies between about 35 Hz and about 60 Hz. 

(13) The apparatus as described in Claim 7 wherein said low-pass filter comprises a standard 
fourth-order low-pass filter at about 50 Hz 

(14) The apparatus as described in Claim 7 further comprising notification means for providing 
notification that a seizure is oncoming in the patient, the notification means being 
communicably connected to said dete rmina tion means 
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